Limiting temperatures and the Equation of State of Nuclear Matter 
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From experimental observations of limiting temperatures in heavy ion collisions we derive T c , 

the critical temperature of infinite nuclear matter. The critical temperature is 16.6 ± 0.86 MeV. 

Theoretical model correlations between T c , the compressibility modulus, K the effective mass, m* 

i i 

and the saturation density, p s , are exploited to derive the quantity (K/rn*) 2 p s 3 . This quantity 
together with calculations employing Skyrme and Gogny interactions indicates a nuclear matter in- 
compressibility in moderately excited nuclei that is in excellent agreement with the value determined 
from Giant Monopole Resonance data. This technique of extraction of K may prove particularly 
useful in investigations of very neutron rich systems using radioactive beams. 

PACS numbers: 24.10.i,25.70.Gh 



Improved knowledge of the nuclear equation of state 
and a coherent picture of the relationship between the 
properties of finite nuclei and bulk nuclear matter re- 
mains a key requirement in both nuclear physics and 
astrophysics. It is key, for example, to understanding 
nuclear structure, heavy ion collisions, supernova explo- 
sions and neutron star properties [[j] ||, gj. Significant ef- 
fort has been devoted to the development of microscopic 
theoretical models which can provide reliable mathemat- 
ical formulations of this equation of state ||, q, 0, [|, [|, 
|| H || II H 0, i@ p. Such calcula- 
tions are usually specified for symmetric nuclear matter, 
a hypothetical system of equal numbers of neutrons and 
(uncharged) protons interacting through nuclear forces. 
Driven by the astrophysical problems and more recent 
laboratory excursions into the region of more exotic nu- 
clei, the dependence of the equation of state on neutron- 
proton asymmetry has also become a subject of signifi- 
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[. In this letter we employ data 
from experimental measurements of caloric curves in nu- 
clear collisions, together with systematic trends and cor- 
relations derived from a number of theoretical investiga- 
tions of nuclear matter, to derive the critical temperature 
and incompressibility of symmetric nuclear matter. The 
techniques employed offer a natural method to extend 
such investigations to more asymmetric systems. 

In a recent paper measurements of nuclear specific 
heats from a large number of experiments were em- 
ployed to construct caloric curves for five different regions 
of nuclear mass [p4[ . Within experimental uncertainties 
each of these caloric curves exhibits a plateau region at 
higher excitation energy, i.e., a "limiting temperature" is 
reached. In Figure 1 these limiting temperatures from 
reference |24|] are presented as a function of mass. As 
previously noted, they are observed to decrease with in- 
creasing mass. This decrease with increasing mass has 
long been predicted as resulting from Coulomb Instabili- 
ties of expanded and heated nuclei p5|, p6| , p7[ p8 , 29, 3(], 



The results employed in reference [[24j were based upon 
temperature determinations derived from double isotope 
yield ratios and from slope measurements of particle spec- 
tra. More recently the TAPS Collaboration has reported 
temperatures determined from a new technique, obser- 
vations of "second chance" bremsstrahlung gamma ray 
emission for a series of reactions which span a wide range 
of massf36], [37| . There are not yet sufficient data of 
this latter type to construct caloric curves for relatively 
narrow mass regions as was done for the previous tem- 
perature data. However, in each case studied with this 
technique the collisions lead to excitation energies which 
are above those identified as the starting points of the 
plateau regions identified in reference) 24 . Thus it is rea- 
sonable to compare the temperatures determined from 
the thermal bremsstrahlung measurements with the ear- 
lier limiting temperature values. As seen in Figure 1, the 
reported second chance gamma temperatures and their 
mass dependence are in excellent agreement with the ear- 
lier results. We take this agreement as an independent 
confirmation of the earlier results and note that the new 
results extend the determination of the mass dependence 
to significantly higher mass. 

A relatively large number of theoretical calculations of 
the critical temperature of semi- infinite nuclear matter 
(nuclear matter with a surface) have been reported in the 
literature g [§ || g |f, || g || @ || . The 
different nuclear interactions employed in these calcula- 
tions lead to large differences in the critical temperatures 
derived. Values from 13 to 24 MeV are reported in ref- 
erences [21 
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limiting temperatures plotted in Figure 1 are well below 
these calculated critical temperatures. This difference 
reflects finite size effects, Coulomb effects and isospin 
asymmetry effects for the finite nuclei studied. A first 
order estimate of the magnitude of these combined ef- 
fects can be made by comparing the volume coefficient of 
the Liquid Drop Model Binding Energy Equation, -16.0 
MeV, which represents the binding energy per nucleon 
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in infinite nuclear matter, to typical nuclear binding en- 
ergies, w 8 MeV/nucleon. This suggests that limiting 
temperatures in nuclei should be « 0.5 times the critical 
temperature of nuclear matter || . Given the wide varia- 
tion in the calculated values of T c it is not surprising that 
large variations result for the absolute values of limiting 
temperatures calculated for finite nuclei. 

Employing a variety of Skyrme type interactions Song 
and Su have previously noted a mass dependent scaling 
of the Coulomb Instability temperatures with the critical 
temperature of nuclear matter (see Figure 6 of reference 
p8| ). These calculations were performed for nuclei along 
the line of beta stability. A similar scaling exists when 
other model interactions are employed. 

Mean values of Tu m /T c for five different masses which 
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27 



the results of different calculations 

H II H HI arc shown in Fig- 
ure 2. The estimated uncertainties are relatively small, 
w 6%. For comparison, the figure also presents ratios of 
Tiim/T c which are expected to result assuming only finite 
size effects as derived from a lattice calculation |58) and 
the ratio of the nuclear binding energy per nucleon along 
the line of beta stability to the bulk binding energy per 
nucleon, 16 MeV. We have employed the mean variation 
of Tn m /T c with A, determined from commonly used mi- 
croscopic theoretical calculations, together with the five 
experimental limiting temperatures reported in reference 
p4| , to extract the critical temperature of nuclear mat- 
ter. In doing so we treat the theoretical variation as if it 
were an experimental uncertainity. Since the various in- 
teractions employed have been "tuned" to other nuclear 
properties, we consider this a reasonable approach. The 
results are presented in Figure 3. Averaging the individ- 
ual results we find 16.6 ± 0.86 MeV. 

It is interesting to ask whether additional Equation 
of State information can be extracted from this result. 
Blaizot et a/.0 have argued that the most effective way 
to extract the incompressibility modulus of nuclear mat- 
ter from experimental data is by comparison with micro- 
scopic calculations. This is usually done by comparison 
of the measured energies for the centroids of the strength 
distribution of Giant Monopole Resonances (GMR) with 
the calculated centroids. The generally accepted best 
current value of K = 231 ± 5 MeV has been determined 
in such a fashion p9| by comparison with the calculated 
centroids using Gogny interactions B. We have adopted 
a similar comparison procedure using the present result 
for T c determination. 

We began by using a relation suggested by the work 
of Kapusta|40| and Lattimer and Swesty |^l| who have 
pointed out that correlations between parameters used 
to describe nuclear matter are such that a relationship 
between the critical temperature, T c , the incompressibil- 
ity, K, the effective mass, m* {— m e ff/m where m e // 
is the nucleon effective mass and m is the nucleon mass) 
and the saturation density, p s , may be written as 



T c = C T (K/m*)2 Ps 



(1) 



where Cy is a constant. Using this relationship, we 
have determined the constant Ct in this equation us- 
ing published theoretical values for T c calculated uti- 
lizing a number of different microscopic interactions 

Results of calculations using interactions with 155 < K < 
384 MeV are depicted in Figure 4. A least squares fit to 
these data suggests a very slight decrease of Ct with in- 
i — - 

creasing (K/m*p p s 3 . Using the present value for T c 
and an iterative technique to establish Ct leads to 



and 



C T = 0.484 ± 0.074 



{K/m*Yp s 3 = 34.2 ± b.UMeV^fm. 



The saturation density, p s , is well established by charge 
density measurements to be 0.16 ± 0.005 fm~ 3 |ll|. The 
standard deviation of the model values of p s from 0.16 
fm~ 3 , calculated with the different interactions, is 3.4%. 
Either of these uncertainties is very small compared to 
other uncertainties in the determination. Therefore T c 
is a measure of K/m*. It is important to recognize 
that K and m* are not independent variables but are 
correlated ||, [13]. For Skyrme effective interactions Cha- 
banet et al. have given expressions for K and m* and 
discussed the resultant correlation between them [Ol . 
Using the relationships discussed in that work m* can 
be written as: 



/ III: ^ I 

= (1 + m p s) 



where 



K-B -Co 



D(l 



la) 



(2) 



(3) 



and a is a parameter which ranges from to 1 and 
controls the density dependence of the interaction. B, 
C, D are parameters directly related to eoo, the energy 
per nucleon in infinite nuclear matter and e f , the Fermi 
energy of infinite nuclear matter. 
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C = -9e c 
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D = T^P k F 



(4) 
(5) 
(G) 



Here kp is the Fermi momentum. 
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It is clear then that for a given K the ratio K/m* in 
equation 1 depends on the choice of a, the parameter 
of the density dependent term. (In the Gogny interac- 
tions of reference |9| this parameter controlling the density 
dependent term is designated a.). As a result, determi- 
nation of K from K/m* is sensitive to the choice of this 
parameter. For example in reference |l3| the relation be- 
tween K and m* is such that small values of a dictate 
lower values of K. Also, for smaller values of a, m* de- 
creases as K increases while for larger tr, m* increases 
with increasing K.(See Fig 2, Reference . 

For comparison to the data we present in Figure 5 a 
i _ i 

plot of K vs (K/m*) ^ p s 3 , obtained using the Gogny 
interactions from reference |^ and various Skyrme inter- 
actions. The dashed lines in the plot show the trend of 
the generalized Skyrme interactions for a — g, 3 and 1, 
as obtained in reference pH} . The solid lines connect re- 
sults for Gogny interactions with a = ^ and |(^] . As seen 
in the figure, a higher value of a leads to a higher appar- 
ent K. It has been pointed out that maintaining K in a 
"reasonable" range of 200-300 MeV requires low values of 
a (or a) ^ |^, [ll| . In particular, the value of K = 231 ± 5 
MeV derived from the GMR data|3{| was obtained by 
comparison of data for the breathing mode energy of 
five different nuclei with energies calculated employing 
the Gogny Dl (a = ±), DlS(a = §) and D250(a = §) 
interactions ||. For three of these nuclei only the D1S in- 
teraction results were used. For the other two a fit to the 
trend in energies calculated from the three interactions 
was employed. 

1 _i 

The value of (K/m*)? p s 3 derived from this work is 
also indicated on the Figure by the vertical line. We 
note that this line intersects the calculated values es- 
sentially at a point where the | Skyrme and | Gogny 
lines intersect. The different slopes of the Skyrme and 
Gogny lines in Figure lead to different uncertainties in 
the K value. Thus employing Skyrme interactions with 
the a = i parameterization fll, K = 232±22MeV. Using 



Gogny interactions with a = 3 |9[ leads to K = 233 ± 39 
MeV. These results for K lead respectively to m* values 
of 0.6741;^ or 0.674±°^. The compressibility modulus 
determined from the critical temperature in this manner 
is then entirely consistent with that determined from the 
GMR measurements. Higher values of a (or a) will lead 
to higher apparent K. Thus for the Skyrme a = | line 
a value of K = 252 would result. For the extension of 
the Gogny a = | line, K = 242 would be obtained). 
The calculated breathing mode energies are apparently 
less sensitive to the value of the parameter of the density 
dependent interaction. 



Extracted by comparison with the same interactions as 
were employed to determine K from observations of the 
Giant Monopole Resonance at low excitation energy, the 
value of K, obtained here from properties of nuclei at 
moderate excitation energies, is found to be in excel- 
lent agreement with that GMR result (39). The precision 
of the GMR measurement is better than that obtained 
from the present determination which incorporates data 
from a number of different experiments. The precision 
for the T c measurement could be improved. However, 
given the relative complexity of the collision dynamics in- 
volved, the breathing mode measurements should remain 
as the standard. Nevertheless, using newly available ra- 
dioactive beams the determination of limiting and criti- 
cal temperatures may play a significant role in providing 
a means to establish the N/Z asymmetry dependence of 
the compressibility modulus and other important nuclear 
properties |42[ . 
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FIG. 1: Limiting Temperatures vs Mass. Limiting tem- 
peratures derived from double isotope yield ratio measure- 
ments are represented by solid triangles. Temperatures de- 
rived from thermal bremsstrahlung measurements are repre- 
sented by open squares. Lines represent limiting temperatures 
calculated using interactions proposed by Gogny (dashed) J3c| ] 
and Furnstahl et al. [B5J (solid). 




FIG. 2: Theoretical variation of the ratio Th m /T c with mass 
along the line of beta stability. The solid line indicates the 
reference value of T c . The short dashed line shows the effect 
of finite size scaling derived from an Ising model [gsj. The line 
with alternating short and long dashes depicts the ratio of 
the nuclear binding energy per nucleon to the bulk binding 
energy per nucleon, 16 MeV. Points with uncertainties are 
derived from the model calculations in references |25[ [ibj |27| , 

H g, |], § §, §,[§. 
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FIG. 3: Derived values of the critical temperature of sym- 
metric nuclear matter. Values derived from data in five dif- 
ferent mass regions are presented. The mean value of 16.6 
MeV is indicated by the horizontal solid line. The range cor- 
responding to ± one standard deviation from this mean value 
is shown by the thin dotted lines. 




FIG. 4: The constant Ct of equation (1), evaluated 
from various microscopic calculations. Ct is plotted against 

(K/m*) 2 p g 3 . Derived values of Ct are indicated by open di- 
amonds. The dotted horizontal line indicates the mean value 
of Ct- The solid line represents the linear least squares fit 
to the derived values. Values of Ct obtained from Skyrme 
and Gogny interactions are further identified by open squares 
placed around the diamonds. 
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FIG. 5: Compressibility modulus, K, as a function of 

j _i 

(K/m*)2p s 3 . The values obtained for Gogny interactions 
of reference |^ are represented by open squares (a = i) and 
open triangles (a = |). Symbols for each set are connected 
by thin solid lines. Results using different Skyrme interac- 
tions are represented by open circles. The other lines repre- 
sent generalized calculations using Skyrme interactions Jl^] 
with a — | (short dashed line), a = ~ (dotted line) and 
a = 1 (alternating dashes and dots). 



